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SUMMARY 

Exponential free-radical growth constants have been measured for ethane- 
carbon monoxide-oxygen mixtures by monitoring the growth of oxygen atom con- 
centration as manifested by CO flame-band emission. Data were obtained over 
the temperature range of 1200 to 1700 K. 

The data were analyzed using an ethane oxidation mechanism involving seven 
elementary reaction steps. Calculated growth constants were close to experi- 
mental values at lower temperatures, up to about 1400 K, but at higher tempera- 
tures computed growth constants were considerably smaller than experiment. In 
attempts to explain these results additional branching reactions were added to 
the mechanism. However, these additional reactions did not appreciably change 
calculated growth constants. 


INTRODUCTION 

Ethane is a simple hydrocarbon and is also an intermediate in the rich 
combustion of other hydrocarbons. Thus, In order to obtain a complete set of 
chemical reactions and rates for modeling hydrocarbon combustion it is desir- 
able to know the pathways by which ethane reacts. 

The shock tube has proved useful for obtaining elementary reaction rates 
from growth constants in the H 2 -CO-O 2 (ref. 1 ) and CH 4 -CO-O 2 (ref. 2 ) systems. 
Consequently an analogous study of the C 2 H 6 -CO-O 2 system was initiated in hopes 
of extracting elementary reaction rates for ethane oxidation. 


Ethane Oxidation Scheme 

When a mixture containing carbon monoxide, oxygen, and a small amount of 
ethane is subjected to a temperature and pressure pulse in a shock tube, small 
concentrations of atoms and free radicals may first be formed by processes such 
as 

CO + O 2 "* CO 2 + 0 
C2H5 + 02"* C 2H5 + HO2 


V\ \ '2 fV- 



and 


followed by 


C 2 H 6 + M -» 2CH 3 + M 


CH 3 + 0 2 ■* CH 3 0 + 0 -» CH 2 0 + 0 + H 

These atoms and free radical concentrations may then grow exponentially via the 
branched chain scheme 


H + 0 2 — ► OH + 0 


OH + CO— ► C0 2 + H 


OH + C 2 H & — ► H 2 0 + C 2 H 5 


H ♦ C 2 H 6 - H 2 ♦ C 2 H 5 


0 + C 2 H 6 — ► OH + C 2 H 5 


C 2 H 5 + °2 C 2 H 4 + H0 2 


10 

C H + M — ► C_H. + H + N 
2 5 2 4 


( 2 ) 

(5) 

(6) 

(7) 

(8) 
(9) 

( 10 ) 


Reactions have been numbered to be compatible with notation used In the study 
of the H 2 -0 2 -C0 system (ref. 1) 

Theory shows (refs. 1 to 5) that In such chain branched systems when 
depletion of reactants Is negligible and temperature and pressure are constant, 
the atom and radical concentrations Increase proportional to exp(Xt) (except 
very early In the reaction). The growth constant X depends on the rate con- 
stants of the elementary chemical reactions and the concentrations of stable 
reactants - In this case CO, 0 2 , and C 2 Hg. 


This report presents experimental growth constants measured behind Inci- 
dent shocks for a range of pressures, temperatures, and gas compositions. 
These experimental growth constants are compared with values computed using 
rate constants from the literature. Calculated growth constants are close to 
experimental values at temperature below 1400 K, but at higher temperatures 
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computed growth constants are considerably smaller than the experimental 
values. Attempts to explain this result using additional chain branching 
reactions have not been successful. 


EXPERIMENTAL ASPECTS 

Growth constants were obtained by measuring the blue CO flame band emis- 
sion behind Incident shocks. The Intensity of this radiation Is proportional 
to the product of carbon monoxide and oxygen atom concentrations (ref. 6), and 
since little CO Is consumed the light measures the Increase of oxygen atom 
concentration with time. Details of the shock tube and associated optical and 
electronic equipment have been described elsewhere (ref. 7). Flame band radia- 
tion centered at 4100 A was Isolated by a high precision ultraviolet Inter- 
ference filter having a 100 A band pass. 

Gas mixtures contained small amounts of C 2 H 6 with varying amounts of 
CO, Og, and C0 2 diluted with argon. (Carbon dioxide was added to ensure 
vibrational relaxation of CO In mixtures with a large CO concentration.) 

Oxygen, carbon dioxide, ethane, carbon monoxide, and argon were high purity 
tank gases and were used without further purification. 

It has been shown (refs. 7 and 8) that boundary-layer effects must be 
considered In analyzing data obtained behind Incident shocks. Conditions 
behind the shocks, In the region of the experimental measurements, were 
obtained from a computer program (ref. 9) which Integrates the equations of 
chemical change for a shocked gas accounting for the effects of boundary-layer 
buildup. The procedure has been described for both turbulent (ref. 7) and 
laminar (ref. 8) boundary layers. The data In this report were obtained with 
laminar boundary layers. 

Exponential growth constants were obtained from plots of logarithm of 
observed light Intensity versus gas time. The relation between gas and labora- 
tory times was obtained from the computer calculations. 


THEORETICAL CONSIDERATIONS 

The analytic solutions of the differential equations describing the Igni- 
tion kinetics In chain-branched systems Involving hydrogen, oxygen, and carbon 
monoxide have been discussed In detail previously (refs. 3 to 5) and the solu- 
tion has been presented for the methane-carbon monoxide-oxygen system (ref. 2 ). 
The solution for the system of reactions presented In the Introduction Is quite 
analogous and will merely be outlined here. 

To obtain the analytic solution, the atom and radical concentrations - H, 
OH, 0, and C 2 H 5 - are assumed to be small In comparison to the concentrations 
of reactants C 2 H 6 , CO, and 0 2 so that the reactant concentrations may be con- 
sidered constant and reactions between chain carriers are unimportant. In 
addition, the temperature and pressure must be nearly constant over the time 
range of the experimental observations. 

Subject to these assumptions, the chemical kinetics are described by a 
system of simultaneous first-order linear differential equations, one for each 
chain carrier. The solution of these equations shows that during the growth 
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period the chain-carrier concentrations grow as exp(xt) (except very early In 
the reaction). During the growth period only reactions (2) to (10) are Impor- 
tant. The growth constant X Is the positive root of a polynomial of a degree 
equal to the number of chain carriers. Thus with H, OH, 0, and C 2 H 5 as chain 
carriers, the following quartlc equation Is obtained: 

-(v 9 + v 10 + x) v 6 

0 -(w 8 + x) 0 

0 v Q - (v 5 + » 6 + X) 

v 10 0 w 5 

where i >2 = k 2 [ 02 l* «5 = ks[C 0 ], = k&[C 2 H&], vj = vg = kgl^Hg], 

vg = kg[ 02 3 , and v-|o = k-jotM]. The rows In this equation derive from the 
differential equations for the rates of formation of C 2 H 5 , 0, OH, and H, 
respectively. The columns derive from the rates of formation or destruction 
of C 2 H 5 , 0, OH, and H In the elementary chemical reactions. 

This equation was used to calculate growth constants corresponding to the 
experimental conditions using rate constants taken from the literature. The 
rate constants and sources are listed In the appendix. Equation (1) was numer- 
ically differentiated to obtain the sensitivities of the calculated rate 
constants to the various reaction rates. 


(v 2 + *•) 


= 0 


( 1 ) 


RESULTS AND DISCUSSION 

The compositions of the three gas mixtures Investigated are given In 
table I. Also shown are the sensitivities, astnx/38.n»i for the various reac- 
tion rates. The range of sensitivities shown Is for the temperature range of 
the data for each composition. Mixture 1 shows sensitivity to the rates of 
reactions 2, 5, 6 , and 8 while mixtures 2 and 3 are sensitive mainly to 
reactions 2 and 8. 

The experimental results are set forth In table II and plotted In 
figures 1 to 3. Because the reactions In the ethane oxidation scheme are all 
blmolecular, the growth constant at a given temperature should be proportional 
to the pressure. Consequently, In figures 1 to 3 the experimental data are 
plotted as the logarithm of x/P versus reciprocal temperature. Also shown 
as solid lines are values computed from equation ( 1 ) using rate constants taken 
from the literature, and set forth In the appendix. 

For mixture 1 calculated growth constants are 25 to 50 percent lower than 
experimental values, with the larger discrepancies at the higher temperatures. 
For mixtures 2 and 3 calculation and experiment are In better agreement, part- 
icularly at temperatures below about 1400 K; but at the highest temperatures 
calculated values are 20 to 50 percent below experiment. 

The agreement between theory and experiment can be Improved by upward 
adjustments In ks and kg. However, the rate of reaction 5 Is quite 
well established, and the Increases required - fivefold or more - cannot be 
justified . 



Other reactions were added to the scheme in hopes of Improving the agree- 
ment between calculated and experimental growth constants. Additional reac- 
tions considered Included 



A. H0 2 + CO -* C0 2 + OH 


B. 


C 2 H 5 + CO -* C 2 H 4 + HCO 
HCO + M -» CO + M + H 

\ 

C. H0 2 + M -» H + 0 2 + M 


D. H0 2 + C 2 H 5 •* H 2 0 2 + C 2 H 5 

j C 2 H5 + 0 2 -* C 2 H50 + 0 
(c 2 H 5 0 + M -» CH 3 CHO + M + H 


The Introduction of another chain carrier - H0 2 or HCO - In schemes A to D 
Increased mathematical complexity by transforming equation ( 1 ) to a qulntlc. 
However, these schemes had negligible affect on calculated growth constants. 


Finally another reaction between atomic oxygen and ethane was added to 
the basic oxidation scheme: 

0 + C 2 H 6 CH 3 O + CH 3 ( 8 a) 

CH 3 + 0 2 -> CH 3 0 + 0 
CH 3 0 + M CH 2 0 + H + M 

Growth constants calculated assuming equal rate for reactions 8 and 8 (a) were 
23 to 24 percent larger for mixture 1 over the experimental temperature range. 
Although this Is an Improvement, It Is not nearly large enough at the higher 
temperatures. If a larger activation energy were assumed for reaction 8 (a) so 
as to Increase the rate of change of growth constant with temperature, the rate 
constants for reaction 8 (a) would be smaller and the Increases In calculated 
growth constants would be less than 24 percent. 

Thus a completely satisfactory theoretical analysis of the experimental 
growth constants has not been achieved. Perhaps the assumptions of the 
linearized analysis - negligible depletion of reactants, negligible heat 
release, negligible radical-radical reactions - are not valid at the higher 
temperatures. 
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APPENDIX 


SPECIFIC REACTION RATE CONSTANTS 

This appendix lists Individual reactions used In calculations and Indi- 
cates values and sources of the assumed rate constants. Blmolecular rate con- 
stants are expressed In cm 3 mole -1 sec -1 . Temperatures are In degrees 
Kelvin and activation temperatures, E/R, are In degrees Kelvin. 

H + 0 2 -* OH + 0 k 2 = 1.38X10 14 exp^^y^) (2) 

This rate constant, from reference 8, was measured In the shock tube used In 
this Investigation. 

OH + CO -> C0 2 + H k 5 = 10 12 exp^y^j (5) 

This rate constant (ref. 1) was determined In the shock tube used In this 
Investigation. 

OH + C 2 H 6 -> H 2 0 + C 2 H 5 k 6 = 8.7X10 13 exp^y 1 ^ (6) 

Rate constant from reference 10. 

H + C 2 H 6 -» H 2 + C 2 H 5 
Rate constant from reference 10. 

0 + C 2 H 6 -» OH + C 2 H 5 
Rate constant from reference 11. 

c 2 h 5 ♦ o 2 -> c 2 h 4 + ho 2 

Rate constant from reference 12. 

C 2 H 5 + M * C 2 H 4 + m + h k io = 10 "exp ^ j (10) 

Rate constant from reference 13. 


, . , 0 14 /- 4881 \ 

k 7 = 1.32x10 exp I — j — I 


(7) 


kg = 1 .1 3x1 0 1 4 expf^y 50 


( 8 ) 


. . c _1 2 / -2445 

kg = 1.5x10 expl — j — 


(9) 


_ i 


mi 
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TABLE I. - MIXTURE COMPOSITION AND GROWTH CONSTANT SENSITIVITIES 


Composition, 
mole percent 

Gas mixture 

] 

1 

2 

3 

C 2 H 6 

0.0215 

0.0190 

0.0204 

CO 

3.04 

10.33 

10.08 

0 2 

9.90 

2.87 

1 .034 

C0 2 

2.74 

5.05 

5.04 

Calculated 




sensitivities, 




reaction 1 




2 

0.19 to 0.11 

0.50 to 0.37 

0.56 to 0.49 

5 

.25 to .26 

.10 to .17 

.06 to .10 

6 

n .09 to .15 

.02 to .03 

.005 to .02 

7 

-.005 to -.001 

-.02 to 1.001 

-.008 to -.001 

8 

.46 to .47 

.41 to .43 

.38 to .40 

9 

-.10 to .02 

-.07 to -.006 

-.015 to -.003 

10 

.12 to .03 

.09 to .011 

.025 to .006 











TABLE II. - EXPERIMENTAL RESULTS 


Temperature, 

K 

Pressure, 

atm 

Growth constant 
X, sec -1 

Mixture 1 

1699 

0.1642 

17.11xl0 3 

1683 

.1421 

17.02 

1677 

.1485 

17.41 

1644 

.1590 

18.48 

1617 

.1584 

16.28 

1605 

.1292 

13.28 

1603 

mmSm: 

12.06 

1580 

Wmmm 

10.01 

1483 

wBuM' 

9.29 

1472 

WeEm* 

6.73 

1409 


7.63 

1390 

.1092 

4.64 

Mixture 2 

1602 


8.82X10 3 

1558 

lKlIlm. 

8.20 

1535 

.1417 

10.52 

1529 

.1510 

9.32 

1498 

.1225 

5.91 

1477 

.1310 

7.15 

1447 

.1229 

5.06 

1438 

.1252 

4.31 

1409 

.1366 

3.79 

1377 

.1188 

3.91 

1359 

.1244 

3.07 

1337 

.1152 

3.45 

1315 

.1170 

3.18 

1287 

.1136 

2.46 

1221 

.1064 

1 .759 

1208 

.1035 

1 .918 









TABLE II. - CONCLUDED 


Temperature, 

Pressure, 

Growth constant 

K 

atm 

\, sec -1 


Mixture 3 


605 

0.1323 

8 . 59 X 10 3 

565 

.1396 

8.25 

536 

.1603 

7.68 

529 

.1278 

4.94 

501 

.1501 

6.03 

463 

.1173 

2.92 

443 

.1272 

4.01 

438 

.1591 

4.91 

432 

.1270 

3.15 

416 

.1406 

3.90 

408 

.1363 

3.07 

388 

.1459 

2.69 

355 

.1157 

1.803 

342 

.1338 

2.12 

321 

.1301 

2.39 

320 

.1201 

2.34 









1700 1600 1500 1400 1300 

Temperature, K 

Figure 1. - Experimental and computed growth constants 
as a function of temperature. Mixture 1. 
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Figure 2. - Experimental and computed growth constants as a 
function of temperature. Mixture 2. 
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Figure 3. - Experimental and computed growth constants 
as a function of temperature. Mixture 3. 
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